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ABSTRACT 


Theoretical evolutionary sequences of models for stars of 

.15 and 30 M have been computed from the zero-age main sequence 

to the end of core helium burning, with the use of the Ledoux 

criterion for convective instability. The sensitivity of the models 

to the most important physical input parameters has been tested. 

During the earliest stages of core helium depletion, the envelope 

rapidly expands into the red-supergiant configuration because the 

hydrogen profile in the intermediate zone remains essentially 

unaltered by local convection. At 15 M , a blue loop on the H-R 

© 

diagram ensues if the initial metals abundance, initial helium 
_ _ 12 

abundance, or C+a reaction rate is sufficiently large, or if 
the 3a reaction rate is sufficiently small. These quantities affect, 
respectively, the opacity (and hence depth) of the base of the outer 
convection zone, the mass of the core, and the thermal properties 
of the core. The blue loop occurs abruptly and fully developed 
when the "critical" value of any of these quantities is exceeded, and 
the effective temperature range and fraction of the lifetime of core 
helium burning during the slow phase of the blue loop vary surpris- 
ingly little. For the recently revised 3a reaction rate, the blue 
loop always forms. At 30 M. , no blue loop occurs for any reasonable 


set of input parameters — unless the temporary convective inter- 
mediate zone (ignored in the present work) extends deeper than doe 
the outer convection zone. For both masses, all passages of the 
Hertzsprung gap occur on the thermal time scale of the envelope. 


I. INTRODUCTION 


The internal structure of stars evolving on the upper main 
sequence is still poorly known (Stothers 1970). Two persistent 
problems are (1) the correct criterion to use for convective in- 
stability in a region that contains a gradient of mean molecular 
weight and (Z) the nature of any resulting convective motions. The 
two possible criteria for convective instability are the Ledoux 
criterion (which includes the gradient of mean molecular weight) 
and the Schwarzschild criterion (which does not). Even were the 
correct criterion known, it is not certain what the final state of 
the convective motions would be. At present, the most popular 
solutions for the internal structure are based on a shrinking 
convective core, in which hydrogen is continually being depleted, 
an intermediate zone of varying hydrogen-helium content, and an 
outer radiative (or mostly radiative) envelope of zero-age 
composition. Although the differences among the various solutions 
have only slight effect on the main sequence phase, the details of 
the further evolution are very sensitive to them. This is part- 
icuhrly true in regard to the direction of evolution (including the 
occurrence of "loops") on the H-R diagram. 
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Few published evolutionary sequences for massive stars have 
been based (implicitly or explicitly) on the Ledoux criterion, and 
have been carried through to the end of core helium burning; those 
available are listed in table 1, where earlier sequences are omitted 
in which the sole opacity source is electron scattering (see Stothers 
and Chin 1968 and references therein). In the present paper, we 
propose to examine in more detail the evolution of massive stars 
based on the Ledoux criterion. 


II. EVOLUTIONARY SEQUENCES BASED ON THE LEDOUX CRITERION 


Theoretical evolutionary sequences for stars with masses of 15. 
and 30 Mohave been calculated by use of an automatic computer pro- 
gram whose basic ingredients are given in the Appendix. Table 2 
lists the following input parameters for each sequence: X^, the 

initial hydrogen abundance by mass; Z , the initial metals abundance 

by mass; the reduced cc -particle width of the 7.12 - MeV level 

'0. 

in ^O; r t the rate of nuclear energy generation from the 3a process; 
a, the ratio of convective mixing length to density scale height; and 
H, the (low-temperature) opacity. The recently revised 2P- rate 
(Austin, Trentelman, and Kashy 1971) will be called here the "new" 
rate, whereas the previously used rate (Clayton 1968) will be called 


the "old" rate. Our basic set of opacities (Appendix) will be referred 



to as the "unmodified" opacities, while the same set multiplied by 
a temperature -dependent increment to the low-temperature (<10000 
°K) values will be referred to as the "modified" opacities. 

Magnetic fields and mass loss have been ignored in our models as 
being probably unimportant (Stothers 1972b_, c_ and references 
therein); rotation has also been ignored although it is undoubtedly of 
some importance. Each sequence has been evolved from the initial 
stage of core hydrogen burning (ZAMS) to the end of core helium 
burning. _ 

Although the general features of evolution in massive stars 
have already been described elsewhere (e. g. , Stothers 1966; 

Stothers and Chin 1968), we shall here summarize the salient 
characteristics of our models because (1) they differ in a few 
significant aspects from our earlier models, (2) they provide a 
framework for our subsequent discussion of changes in the input 
parameters, and (3) they are relevant to our discussion of the 
relative amounts of time that a massive star spends burning core 
helium in various parts of the H-R diagram. Three representative 
evolutionary tracks for each mass, based on three different initial 
chemical compositions, are shown in figure 1. The interior dis- 
tributions of hydrogen at several stages of evolution for two of 
these initial chemical compositions are shown in figure 2 for 15 M . 



Beginning at the ZAMS, the main phase of core hydrogen 

burning is characterized by an extensive convective core, which 

depletes its hydrogen, recedes monotonically in time, and leaves 

behind a convectively stable intermediate zone of varying hydrogen- 

helium content. The slowly expanding outer envelope is also 

convectively stable, except for a small region just above the zero-age 

core boundary in 30 where convection develops during the late 

1 

stages of core hydrogen burning. 

■^If the adopted opacity source is pure Thomson 
scattering, convection develops just outside the zero-age 
core for masses as low as 10 (Schwarzs child and Harm 

1958). Furthermore, if convective overshooting into the 
immediately underlying radiative layers is important, a 
small semiconvective zone will develop there (Sakashita 
and Hayashi 1961; Stothers 1963). ____^ 

When the central hydrogen content drops to X c = 0. 030 ± 0. 005 
(see also Stothers 1964), the star undergoes an overall contraction; 
the hydro pen immediately surrounding the convective core heats 
up and a hydrogen-burning shell is formed.. 

After the complete exhaustion of hydrogen at the center, the 
core undergoes a rapid contraction, while the envelope re-expands 


vigorously. Convection at the center eventually vanishes during 
this phase if the stellar mass is less than~30 2 (see also 

~Z " ’ ' 

More exactly, if the core mass is less than 

^*0 The corresponding mass for an ideal core 

characterized by a pure helium composition, pure 

Thomson scattering opacity, and purely homologous 

contraction is 8 M q (Tanaka and Sakashita 1964). 

i 

Stothers 1964). Ignition of central helium causes the redevelopment 
of central convection, and helium burning finally stabilizes the 
contracting core just before, or during, the crossing of the 
Hertzsprung gap. However, the envelope keeps on expanding (on 
a thermal time scale) until the tip of the red- supergiant branch is 
reached. By this time, some 3 percent of the core helium is 
consumed. 

The region immediately below the zero-age core boundary 
becomes convectively unstable during the crossing of the Hertzsprung 
gap. Any temporary convective modifications of this region (ignored 
in our models but included by Stothers 1966) are soon obliterated by 
the deep convective envelope which extends down from the surface 
“certainly this happens at 15 M although possibly not at 30 M 



(see Chiosi and Summa 1970; Robertson 1972). The maximum 


depth achieved by the convective envelope occurs when the star is 
at or near the tip of the red-supergiant branch. There the surface 
hydrogen abundance is reduced by 1 to 10 percent (depending on 

the mass and initial chemical composition ) because convection 

/ 

penetrates to the helium-rich intermediate layers; simultaneously, 
12 14 

the surface C/ N ratio is also reduced (see Iben 1964) although 
we have not determined the amount here. 

Subsequent evolution down the red-supergiant branch takes 
place slowly on a nuclear time scale. The descent is halted 
prematurely if the envelope becomes again thermally unstable and 
contracts rapidly to the blue-supergiant configuration; thereupon, 
depletion of core helium proceeds within a thermally stable blue 
supergiant. Sometimes a small S-bend, or stillstand, occurs at 
the blue tip of the evolutionary loop, before the envelope re -expands 
vigorously during the last stages of core helium burning, ultimately 
to cross again the Hertzsprung gap on a thermal time scale (usually 
after the exhaustion of central helium). However, if the star's 
envelope on the red- supergiant branch remains thermally stable 
throughout the phase of core helium burning, the descent of the red 
branch continues to fainter luminosities, down to a precarious point 
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where the star may "hover" for a time before making a final 

re-ascent. Envelope convection at this minimum luminosity is 

found to cover a much smaller mass fraction at 15 M than at 30 M 

0 © 

2 _ 2 

-^ or Q “ 0 an< f 9 - 1> the final chemical composition of the 

CL Ot 

helium-exhausted convective core is always pure carbon and pure 

2 

oxygen, respectively. For ^ = 0. 1, the final carbon abundance is 
X = 0. 15 ± 0. 05 for the new 3a rate (or X = 0. 35 ± 0. 05 for the 
old 3a rate). 

The last seven columns of table 2 give the following data: 

l°g logarithm of the effective temperature at the tip of the 

blue loop during core helium burning; T , lifetime of core hydrogen 

burning (until X c =0); ratio of the lifetime of core helium 

burning (until Y = 0) to that of core hydrogen burning; T (1)/t , 

c b He 

T b^ 2 ^ T He ’ fraction of time spent in the first (1) and second (2) 
blue - super giant configurations (log ^ 4. 1) during core helium 
burning; t^(1)/t^ g , T y(2)/T ^ > fraction of time spent in the first 
(1) and second (2) yellow-supergiant configurations (3. 6 < log T 

e 

^4. 1) during core helium burning. The last significant figures 
given for me blue and yellow lifetimes are rather uncertain 
because of the uncertainty of the effective temperature boundaries. 
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Generally speaking, our evolutionary sequences are in good 
qualitative agreement with previously published sequences based 
on the Ledoux criterion (table 1). Our specifically new results, 
and some differences of interpretation, are discussed below. 

III. INVESTIGATION OF 15 M 

0 

Our most extensive set of model sequences, designed to test 
the effect of changing various physical parameters, has been 
carried out for 15 M^. Properties of the sequences, other than 
the blue loops, which are discussed below, will be analyzed next 
in terms of these physical parameters. 


a) Nuclear Reaction Rates 


Hydrogen burning via the CN cycle completely outweighs 

the pp chain, which has been ignored in our models. However, 

14 

possible errors in the N + p reaction rate and, more significantly, 

14 

in the assumed equilibrium abundance of N have been assessed 

1 2 

by considering the two cases X - — Z and X - — Z , 

UIN j 6 UlN U j G 


which are introduced as the N abundance in the reaction rate. 


Virtually no difference in results is found between our two cases 
(table 2) except for a slight lengthening of lifetimes and a slight 
brightening of the tip of the red -supergiant branch due to enhanced 



shell burning and to an enlarged core mass, when N is more 
abundant. 

Nitrogen burning via a capture has been ignored as an energy- 

source in all of our models since the revised reaction rate for 

14 3 

N + tt implies an unfavorable competition with 3a. However, the 

3 14 

Even with the old rate, N + a scarcely affected 

the model sequences {Iben 1966a, b; Simpson 1971). 


revised reaction rate for 3a causes helium to burn at a hotter tem- 
perature than before, and so produces a somewhat brighter tip of 
the red- supergiant branch and a deeper convective envelope; the 

elevated temperature also means that carbon is depleted faster 

• 12 „ . , .12 

via C + a . At the present time, C + a can be considered the 

only important reaction whose rate is still uncertain. Model 

Q Z 

sequences based on the choices of a =0, 0. 1, and 1 for this 

a 

Z 

reaction show than an increase in q (by augmenting the formation 
of oxygen as well as the size of the convective core) prolongs the 
lifetime of helium burning. 

b) Convective Envelope Parameters 


Our sequences are found to be not greatly affected by changes 



in the physical parameters a and * pertaining to layers above 
the base of the outer convection zone, because the location of the 
base is itself not much affected. This result has been obtained 
earlier (Stothers 1972a), and is due to the fact that, in massive 
stars, the superadiabatic region near the surface covers only a 
small fraction of the total mass. However, the surface temperature 
can be adjusted to agree with observations by selecting an appro- 
priate value of CL or of low-temperature opacity (Henyey, Vardya, 
and Bodenheimer 1965; Stothers 1972a), although no physical 
significance should be attached to such arbitrary selections, which 
merely mask the uncertainties in the mixing -length theory. 

c) Initial Chemical Composition and Opacity 

The initial hydrogen (and helium) content of the star influences 
its structure primarily through the mean molecular weight (and, to 
a lesser extent, through the opacity). As can be shown by simple 
similarity arguments, an increase in reduces the luminosity 
and mass fractions of composition boundaries (including the base 
of the outer convective envelope), and therefore increases the life- 
time of the various phases. However, as shown by figure 3, the 
mass separations between composition boundaries are nearly 
independent of X , as is the ratio of the hydrogen-shell luminosity 
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to the total luminosity (at a given stage of evolution). But the 
fraction of the total stellar lifetime spent burning core helium is 
slightly less for large values of X g because less oxygen is formed 
in smaller cores. 

The changes induced by an increase in Z are quite similar 

to those induced by an increase in X , with the main exception that 

e 

the convective envelope penetrates to a disproportionately deeper 

level (see figure 3). Such a change can be induced also by an 

increase in the opacity coefficient alone, but this is equivalent to 

an increase in Z since the latter quantity influences the structure 
e 

of the star primarily through the opacity? The depth of the con- 
vective envelope is determined, in part, by the opacity near its base, 

6 

where the temperature (for a deep envelope) is normally ~10 °K; 

if the opacity is large at this temperature, the radiative temperature 
gradient will remain steep to a greater depth than otherwise. 

d) The Blue Loops 

The chief feature which differentiates the various sequences 
on the H-R diagram is the appearance, or nonappearance, of a 
long blue loop during core helium burning. The precise cause of 
the loops is still obscure. The "mirror principle" (Hayashi, Hoshi, 
and Sugimoto 1962; Hofmeister, Kippenhahn, and Weigert 1964), 



which is a phenomenological description (not always valid) of what 
happens, has been analyzed in a preliminary way by Stein (1966). 

The roles of X^, of core mass, of opacity, and of hydgrogen 
profile in and above the shell have been recognized in our previous 
study (Stothers and Chin 1968), while Hofmeister (1967) has 
emphasized the importance of Z . Further studies (Fricke, Stobie, 
and Strittmatter 1971; Lauterborn, Refsdal and Roth 1971; Barbaro, 
Chiosi, and Nobili 1971; Dallaporta 1971; Fricke and Strittmatter 
1972) have confirmed and greatly amplified the earlier work, and 
it now seems that such additional features as numerical accuracy 
in the computations (Paczynski 1970b) and axial rotation (Kippenhahn, 
Meyer-Hofmeister, and Thomas 1970; Meyer-Hofmeister 1972) are 
also important in determining the loops. j 

Whether or not a blue loop occurs seems to depend on the 
existence of multiple solutions of the static stellar configuration 
for a given run of chemical composition (Paczynski 1970b; 

Lauterborn et al. 1971; Kozlowski 1971). 4 It has further been 

4 

Such a multiplicity of solutions was already 
found among (purely) red- supergiant models for more 
highly evolved states (Stothers and Chin 1969). 



shown that initiation of the loop for a massive star takes place on 
the thermal time scale of the envelope. We have found that the 
closer a sequence is to a "critical" sequence (one which barely 
develops a blue loop), the more likely to occur, and the more 
pronounced, is a thermal fluctuation of the envelope just around 
the luminosity minimum before the star goes either over to the 
blue tip or back up to the red tip. To investigate the possibility 
of thermal instability of the hydrogen-burning shell in connection 
with the loop phenomenon, we have recalculated parts of our 
sequence 15 -P using very small time steps, but have found that 
the shell always remains thermally. stable (as expected from our 
analogous study of the earlier stages of evolution [Stothers and 
Chin 1972]). However, the tendency toward thermal instability in 
the shell is greatest during the initial phase of the loop. Since the 
envelope contracts very rapidly at this time, the shell heats up 
faster than it is cooled by its expansion; therefore it temporarily 
spreads over a greater mass fraction, and its luminosity rises. 

This relieves some of the weight off the core, so that the central 
pressure, and hence the central temperature and the convective -core 
boundary, initially diminish a little. Thus, there is feedback from 
the envelope to the core, and the "isolated star" approximation for 
the core structure is evidently only valid for quasi-static stages of 


evolution. 
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Although we are unable to give a precise physical explanation 
for the occurrence of the blue loops, we can, at least, point out 
four factors which lead to their development: (1) a small value of 

X , (2) a large value of Z {or of interior opacity), (3) a large value 

G C 

12 

of the C + CL reaction rate, and (4) a small value of the 3^ reaction 
rate. Moreover, a blue loop, when it occurs, does so abruptly for 
a small change in one of the four parameters, and always appears 
fully developed on the H-R diagram. 

Interestingly enough, core mass alone (as fixed by X ) can 

e 

determine whether or not a blue loop occurs. Our sequences for 

15 M with Z =0. 021, 0 = 0. 1, and the old 3& rate, have a loop 

O e a ^ 

only if X ^ 0. 64. But, regardless of X , they all contain the same 
e e 

hydrogen gradient in and above the burning shell (at a given stage 

of evolution), and, in all cases, the shell never makes contact with 

the hydrogen discontinuity that is formed by the outer convection 

zone (see figures 2 and 3). The effectiveness of core mass can 

also be seen in our previously published sequence for 15 M 

O 

(Stothers and Chin 1968), which initiated a loop at the same core 
mass (q .^0. 32) as do our present sequences — even though the 
hydrogen profile in the shell was very different. 

In contrast to the sequences with low Z , our sequences for 



15 with = 0. 044 have relatively deep convective envelopes. 

Before the end of core helium burning, the hydrogen-burning shell 

reaches the hydrogen discontinuity in all cases (see figures 2 and 3), 

and this produces a long blue loop (see also Stothers and Chin 1968; 

Robertson 1972). It appears likely that a difference in Z can also 

e 

explain why the sequence for 20 calculated by Chiosi and Summa 

has a blue loop whereas Robertson’s sequence, for very nearly the 

same value of X , does not (see table 1). Paczynski's adopted 

metals abundance, Z =0. 03, used in his sequence for 15 M , is 

e 0 

"marginal" in the present context, and this may explain why, in 
two preliminary trial sequences, he obtained a blue loop in one case 
and no blue loop in the second case (significantly, the loop occurred 
at almost the end of core helium burning). Our present sequence 

f 

15-U spends less time in the blue loop than does Robertson's 
analogous sequence, because his adopted opacities were larger, 
resulting in a deeper convective envelope at the tip of the red-super - 
giant branch and hence an earlier conjunction of the hydrogen- 
burning shell with the hydrogen discontinuity. The variation between 
Hofmeister's two sequences for 9 M^is probally due to the difference 
in X , although the difference in Z may be important if the convective 
envelope in her sequence with large Z is sufficiently deep. 



2 

Finally, if 0 =1 and the old 3a rate are adopted, our sequences 

a 

for 15 M with Z =0. 0Z1 exhibit a blue loop whenever X ^ 0. 70 
© e e 

(to be compared with the much smaller limit on X for the corre- 

2 

sponding sequences with 0^ = 0. 1). But, if the new 3a rate is 

adopted, all the sequences for 15 M Q have a blue loop, regardless 

of the adopted values for the other parameters. Since the changes 

12 

in the rates of the C + a and 3a reactions do not alter the mass 
of the core or the hydrogen gradient in the envelope, the changes 
of the thermal properties (alone) of the core are, clearly, the 
determining factor. However, in the case of the new 3a rate, the 
slightly deeper location of the hydrogen discontinuity does, for 
some sequences, assist in promoting the blue loop. 

IV. INVESTIGATION OF 30 M q 

All of our evolutionary sequences of 30 remain on the red- 
supergiant branch during core helium burning. There are two major 
reasons for this. First, the hydrogen gradient, dX/dM (r), in the 
intermediate zone is very shallow at high mass since the non- 
dimensional gradient dX/dq is nearly independent of mass. Second, 
the outer convection zone does not penetrate very deeply into the 
star (figure 4), in contradiction to our earlier prediction (Stothers 
and Chin 1968). Although it is true that convection extends down 


well into the intermediate zone just as it does for lower masses, ^ 

5_ 

For some reason, the outer convection zone did not 
reach the intermediate zone in Simpson's (1971) sequence 
for 30 (see his figure 5a). 

this zone is displaced farther outward in mass fraction at higher 

masses, and, because of the generally shorter time scale, the hydro- 

gen-burning shell does not advance as far as it does at lower masses, 

with the result that it fails to reach the hydrogen discontinuity. An 

earlier sequence calculated by Simon and Stothers (1970) for (X , Z ) 

e e 

= (0. 33, 0.03) did have a blue loop (without the aid of a suddenly 

encountered hydrogen discontinuity), so that the "critical" hydrogen 

abundance for 30 M probably lies somewhere between X = 0. 33 and 

O e 

0.55. The present sequences for 30 are subject to far less un- 
certainty about the opacity at temperatures near 10^ °K than are the 
sequences at lower mass, due to the predominance of electron 
scattering at high mass. To check this point, we have computed a 
special sequence for 30 with all input parameters favorable 

to the development of a blue loop: X = 0. 602, Z = 0. 044, 0 2 = 1 

e e a 

the new 3a rate, and the absorption part of the opacity multiplied by a 
factor max f 1, 3 - | 6 - log T|] (inspired by the opacity calculations 
of Carson, Mayers, and Stibbs 1968). The sequence remained red. 
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Nevertheless, our neglect (and Simpson’s neglect) of semi- 
convective and convective mixing in the intermediate zone may be 
serious for masses heavier than 15 (§11). It is possible that a 

convective intermediate zone might extend to a greater depth than 
that reached by the outer convection zone, and hence induce a blue 
phase when the hydrogen-burning shell reaches the hydrogen discon- 
tinuity. Further studies are necessary to decide on this question. 

Bisnovatyi -Kogan and Nadezhin (1972) obtained results for 
30 that are qualitatively different from Simpson's and our results. 
Since these authors were unable to obtain convergent, static solutions 
for the atomspheres of yellow-supergiant models and also found a 
large density inversion near the surface of their last static models, 
they postulated a steadily outflowing atmosphere as the proper 
boundary condition for the next stages, and derived a mass loss 
rate of 0.5 per year. We object to their solution on three 
counts. First of all, our probably more accurate treatment of the 
static atmospheres (see Appendix) always led to consistent, con- 
vergent solutions ; more importantly, Osmer (1972) also obtained . 
convergent solutions for much more detailed model atmospheres. 
(Both Osmer's and our models contain density inversions, but 
so do models for main-sequence stars.) Second, it is 
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observationally (Osmer 1972) and theoretically (Wentzel 1970) 
likely that a density inversion will lead to turbulent motions that 
adjust the atmospheric structure appropriately rather than to drastic 
mass loss. Third, and finally, the observed masses, mass-loss 
rates, and even existence of red supergiants of high mass are in 
conflict with the idea of heavy mass loss from yellow -supergiant 
progenitors (Stothers 1972b), quite apart from the problem of how 
a massive main-sequence star could have survived its red and 
yellow-pre-main-sequence stages. 

V. CONCLUSION 

Adoption of the Ledoux criterion for convective instability in 

calculating the evolution of a massive star causes the development 

of a shallow hydrogen gradient in the intermediate zone, which makes 

the envelope expand rapidly into the red-supergiant configuration 

during the early stages of core helium depletion. A blue loop on 

12 

the H-R diagram ensues if Z , Y , or the C + a reaction rate 

e e 

is sufficiently large, or if the 3a reaction rate is sufficiently 
small. It appears that these quantities are related, respectively, 
to the influence of opacity on the location of the base of the outer 
convection zone, to the mass of the core, and to the thermal 
properties of the core. The blue loops, however, are not 
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particularly sensitive to the choices of X a. or low-temperature 
opacity. Our main results are summarized in table 3. Whether 
the star crosses the Hertzsprung gap once or three times during 
core helium burning, all of these passages occur on the Kelvin 
time scale of the envelope (approximately). Before each passage, 
a slight pause or fluctuation of the luminosity always develops. 
Exhaustion of fuel at the center, and the consequent switchover to 
a purely gravitational source of energy in the core, is the reason 
for the luminosity perturbation before the first and third passages. 
Thermal instability of the envelope, and the consequent reduction 
of pressure on the core, accounts for the luminosity perturbation 
before the second passage. A general rule concerning the evolu- 
tionary tracks on the H-R diagram is that the fine details become 

more obvious when X or Z is large. 

e e 

Unfortunately, the blue-to-red ratio for supergiants of~15 M 
cannot be predicted with very much precision until accurate know- 
ledge is obtained of 0^ and of the stellar opacity at temperatures 

6 2 

near 10 °K. At present, a rather small value of 0 (~0. 1) is 

a 

expected to be correct (Barker 1971; Vidal, Shaviv, and Kozlovski 
1971), but the stellar opacity is a more perplexing problem 
(especially near 10^ °K) since the relevant atomic physics as 



well as Z and the constituents of Z are not yet well enough known 
e e 

(Carson et al. 1968). Our present models for 15 indicate, 

tentatively, that a blue loop will occur for any reasonable initial 

chemical composition (providing that the new 3a rate is correct but 

2 

regardless of the value of 0 ) and also that j / t t , = 0. 5 ± 0. 2. 

& a b He 

A more precise statement can be made about the maximum 
effective temperature which is attained if the sequence has a blue 
loop. We find that this quantity is relatively invariant for our 
sequences at 15 M^, which agree well, in this respect, with the 
two sequences obtained earlier by Stothers and Chin (1968) and by 

Robertson (1972). ^ Thus, helium burning in blue supergiants of 



6 

However, the maximum effective temperature 
attained is slightly cooler than average for the earlier 
sequence calculated by Stothers and Chin because of their 
specification of no chemical homogenization at all in the 
intermediate zone and their adoption of exaggerated low- 
temperature opacities (given by Kramers law). The opacity 
effect alone can be seen explicity in our present work by 
comparing sequences 15-P and 15-Q. On the other hand, 
Robertson's sequence becomes slightly hotter than average 
because of extensive chemical homogenization in the 
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intermediate zone (due to a deep convective envelope on the 
red-supergiant branch). 

~15 M q (based on the Ledoux criterion) is expected to occur within 

an effective temperature range from log T « 4, 1 to log T - 

c c 

4. 22 ± 0. 03, since stabilization of the rapidly contracting envelope 

7 

always first occurs when log ^ 4. 1. 

7 2 

The effect of increasing the parameters 0 or Z , 

cc ^ 

or of decreasing the 3a rate, is to increase only moderately 
t / t and log T at the tip of the blue loop. The corre- . 
sponding effect of increasing of X^.^. is insignificant. 

At 30 M q , adoption of the Ledoux criterion does not yield a 
blue loop for any reasonable initial chemical composition. Unless 
convection in the intermediate zone (which we have ignored) pene- 
trates to a deeper level than does the outer convection zone, the 
loop is unlikely to occur at this mass. 



APPENDIX 


THE STELLAR -EVOLUTION PROGRAM 

The stellar -evolution program developed for our work has 
been based on the standard method of numerical relaxation; the 
particular version of the method employed is that of Kippenhahn, 
Weigert, and Hofmeister (1967), with only minor modifications. 
We shall present here the details of the input physics that we have 
adopted because such information is useful for comparisons with 
other work, for assessing present limitations of the physics, and 
for justifying our various approximations^ The same program ha: 
already been used for the models in two recent papers of ours 
(Stothers 1972b; Stothers and Chin 1972). f 

a) Thermodynamic State 

The thermodynamic state of the interior (in massive stars) 
is represented accurately by a mixture of nondegenerate gas and 
radiation. The state of the gas is computed for the following 
elements: hydrogen and helium in the stages H^, H, H + , He, He 
and He ; a hypothetical metal, with low ionization potential, in 
the stages M and M + ; and all other metals in the fully ionized 
state or the neutral state (depending, for simplicity, on whether 



helium is or is not completely ionized). The hypothetical metal 

is assumed to have the following atomic properties: 1=7.5 eV, 

B(1)/B(0) = 1.25, A = 40, and = 0.08 Z g (Demarque I960; Iben 

1963). The average atomic weight of all the other components of 

Z is taken to be A = 16 (Traving 1966). 
e 

b) Opacity 

\ 

Opacities calculated by Cox and Stewart (1965), including 
line absorption, have been adopted in the form given by Stothers 
and Simon (1970) as an improvement of Christy's (1966) formula, 
for temperatures higher than 4000 °K. H opacities as given by 
Kippenhahn, Temesvary, and Biermann (1958) have been adopted 
for temperatures below 3000 °K. A temperature -weighted mean 
between the two sets of opacity has been employed for the inter- 
vening temperature range. Compton scattering by free electrons, 

in the form given by Deinzer and Salpeter (1965), has been adopted 

7 

for temperatures higher than 8/10 °K. 

Corrections to this prescription for opacity are undoubtedly 
needed, but the direction and amount of the total correction are 
still very uncertain. For example, compare the cited references 
in regard to the following quantities for opacity at high temperatures 



the atomic model adopted (Carson, Mayers, and Stibbs 1968; 
Nikiforov and Uvarov 1970); the relative abundances of minor 
elements (Carson_et_al. 1968; Ezer and Cameron 1971); auto- 
ionizing states (Watson 1969a; Carson 1972); and electron corre- 
lations (Watson 1969b; Diesendorf 1970). Fortunately, in our 
models of massive stars, the densities are too small for electron 
correlations to be important. At low temperatures, however, we 
have omitted H^, H^O, CO, and CN, which are probably the most 
important molecular sources of opacity; but, even so, for high 
masses these sources can influence the interior structure of the 

stellar envelope very little (§ Illb; Auman and Bodenheimer 1967). 

12 

Our only important omission seems to be absorption due to C 
16 

and O in the dehydrogenized core, but the resulting formal 
error in opacity is ^ 10 percent (cf. Cox and Stewart 1965). 

c) Nuclear Reaction Rates 

The following nuclear processes have been taken into account: 
4 -» 4 He (CN cycle), 3 4 He - 12 C, and 12 C + 4 He 16 0. 

Energy-generation rates for the equilibrium reactions have been 
taken from Clayton ( 1 9 68 ) . For hydrogen burning, we have assumed 
that (Traving 1966) while ignoring the contribution of 

O to the full CNO bi-cycle. The recently revised 3ct rate (Austin, 



Trentelman, and Kashy 1971) is smaller than the rate quoted by 

Clayton by a factor exp{-l. 38/ T ), and has been used in our 

1 2 

latest model sequences. The rate for C + Ct is so uncertain 

2 

(Barker 1971) that we have left as a free parameter 0 , the 

reduced ct -particle width of the 7. 12 - MeV level in O. The 
14 

revised rate for N + a (Couch et ah 1972) indicates that nitrogen 

will be depleted only after the onset of the 3a reactions, and so 
14 

N + a has been neglected as an energy-generation source. 
Electron-screening corrections to all the rates have been ignored 
as being unimportant because of the low densities encountered in 
massive stars. - 

d) Convection 

I 

Convection in the deep interior has been treated as fully 
adiabatic. In the outer envelope, nonadiabatic convection has been 
calculated according to the prescription given by Iben (1963, 1965), 
who adopted Bohm-Vitense 1 s (1958) version of the mixing -length 
theory. Both the low -temperature opacity and the ratio of con- 
vective mixing length to density scale height have been left as 
adjustable parameters. 

Convective overshooting at the inner boundary of the convec- 
tive envelope has been calculated in an approximate way, which 
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mimics a very thin semiconvective zone (Ledoux 1947; Chiosi and 
Summa 1970). At the outer boundary of the helium -burning 
convective core, overshooting has been neglected because of its 
already demonstrated insignificance in massive stars (Robertson 
1972). All boundaries between radiative and convective zones have 
been fixed by requiring equality of the radiative and adiabatic 
temperature gradients on the convective sides of the boundaries 
(Harm and Schwarzschild 1955). 

e) Semiconvection 

Semiconvection has been ignored in our models even though 
it is occasionally implied by the Ledoux criterion for convective 
instability in the intermediate zone of varying hydrogen-helium 
content. Its presence at the inner boundary of the convective 
envelope (above) and at the outer boundary of the helium -burning 
convective core (Robertson 1972) has not been explicitly taken into 
account in our models, since, for massive stars, these semicon- 
vective boundary layers are very thin. 

f) Surface Boundary Conditions 

The uppermost layers of the atmospheres of our models 
have been assumed to be in radiative equilibrium. Eddington's 



(1926) approximate solution for a radiative atmosphere has been 
used for the local temperature structure, and an accurate inte- 
gration of the equation of hydrostatic equilibrium has been performed 
for the pressure structure. The inward atmospheric integration is 
terminated either at the layer where the local temperature is equal 
to the effective temperature of the star or at the layer where con- 
vection begins (whichever comes first). The inward envelope 
integration is then begun as a smooth continuation from this point. 

It is terminated at mass fraction q = 0. 90, where the values of the 
physical variables are taken for use in the relaxation part of the 


program. 
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TABLE 1 

PREVIOUS EVOLUTIONARY SEQUENCES OF MODELS 
FOR MASSIVE STARS BASED ON THE 
LEDOUX CRITERION 


M/M 

O 

X 

e 

Z 

e 

9 2 

a 

e 3a 

tog T e 
(tip) 

T b 

/t TT 

He 

Author 

V 

9 

0.739 

0. 021 

0. 100 

old 

• • • • 

0. 0 

Hofmeister (1967) 

9 + 

0. 602 

0. 044 

0. 100 

old 

4. 03 

0. 5 

Hofmeister (1967) 

10 

0. 700 

0. 030 

0. 085 

old 

4. 11 

0. 5 

Paczynski (1970a, b) 

15 

0. 700 

0. 030 

0. 085 

old 

• ' • • ■ 

0. 0 

Paczynski (1970a, b) 

15 

0. 700 

0. 030 

0. 100 

old 

4. 17 

0. 6 

Stothers and Chin (1968) 

15 

0. 602 

0. 044 

0. 100 

old 

4. 25 

0. 6 

Robertson (1972) 

20 

0. 617 

0. 020 

0. 100 

old 

.... 

0. 0* 

Robertson (1972) , 

20 

0. 602 

0. 044 

0. 060 

old 

4. 27 

0. 8 

Chiosi and Summa (1970] 

30 

0. 750 

0. 020 

0. 085 

old 

• • • • 

0. 0 

Simpson (1971) 


* Also calculated by Kippenhahn, Meyer -Hofmeister, and Thomas (1970) and by 
Meyer-Hofmeister (1972), with and without axial rotation. 

^ Also calculated by Lauterborn, Refsdal, and Roth (1971), by Fricke, Stobie, and 
Strittmatter (1971), and by Fricke and Strittmatter (1972). 

+ 

This sequence was terminated when Y c = 0. 447, but no blue loop had yet 

,/ 

II 


developed. 


NEW EVOLUTIONARY SEQUENCES OF MODELS FOR STARS OF 
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TABLE 3 

INPUT PARAMETERS FOR WHICH 

«A» 

A BLUE LOOP DEVELOPS"' 


M/M 

O 

X 

e 

Z 

e 

9 2 
a 

€ o 

3a 

15 

^ 0. 62 

0. 021 

o 

• 

o 

old 

15 

< 0. 64 

0. 021 

0. 1 

old 

15 

< 0. 70 

0. 021 

1. o 

old 

15 

all 

0. 044 

0. 1 

old 

15 

all 

all 

all 

new 

30 

none 

none 

none 

new or old 

* 'i' 

An 

entry "all" or "none" 

refers to the following 

ranges: 


X 

e 

= 0. 602 - 0. 739, Z 

e 

= 0. 021 - 0. 044, and 0 

2 = 0. 0 - 1. 0. 
a 




FIGURE CAPTIONS 


Fig. 1. — Theoretical evolutionary tracks in the H-R 
diagram for stars of 15 and 30 M^. Tracks are labeled with the 
sequence numbers (see table 2), indicating the differences in 
initial chemical composition. 

Fig. 2. — Hydrogen profiles for stellar models of 15 

with Z =0. 021 (solid lines, sequence 15-0) and Z =0. 044 
e e 

( dashed lines, sequence 15-F). Lettering index: a, end of core 
hydrogen burning; b, tip of the red-supergiant branch; and c_, 
end of core helium burning. 

Fig. 3. — Mass fractions of important composition 

boundaries as a function of X , for the sequences at 15 M with 

e i O 

Z =0. 021 (solid lines) and Z = 0. 044 (dashed lines). Lettering 
e e 

index: a, maximum extent of the hydrogen-burning (ZAMS) 
convective core; b, innermost extent of the outer convective 
envelope; c, hydrogen shell peak at the end of core helium 
burning; d, hydrogen shell peak at the tip of the r ed- supergiant 
branch (during the early stages of core helium burning); and e, 
maximum extent of the helium -burning convective core. 

Fig. 4. — Mass fractions of important composition bound- 
aries as a function of X , for the sequences at 30 M . Notation 


is the same as in figure 3, 
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